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We introduce a dynamical fermion algorithm which is based on the hybrid Monte Carlo (HMC) algorithm, but 
without pseudofermions. The molecular dynamics steps in HMC are retained except the derivatives with respect 
to the gauge fields are calculated with the Z2 noise. The determinant ratios are estimated with the Pade - Z-i 
method. Finally, we use the Kennedy-Kuti linear accept/reject method for the Monte Carlo step which is shown 
to respect detailed balance. We comment on the comparison of this algorithm with the pseudofermion algorithm. 



1. Introduction 

The partition function of the gauge field theory 
with dynamical fermions in the Euclidean path- 
integral formalism has the following generic form 

Z = y" \dU\detM\U\e-^'=' (1) 

where Sq is the pure gauge part of the action and 
U is the link variable. The Grassmann numbers 
representing the fermions are formally intergrated 
out to give a fermion determinant. Since the de- 
terminant is very non-local, the Monte Carlo al- 
gorithms with dynamical fermions is known to 
be very time-consuming. The best algorithm so 
far is the Hybrid Monte Carlo [Q] which intro- 
duces a Monte Carlo step to the hybrid (molec- 
ular/hybrid) algorithm to remove the discretiza- 
tion error from the molecular dynamics step. This 
algorithm avoids calculating the determinant di- 
rectly by introducing pseudofermion variables. 
The partition then takes the form 

J[dU][dx*][dx]e~^^''''^'''''^"'' (2) 

where x is the pseudofermion field variable. 

In this paper, we shall present an algorithm 
which estimates the determinant ratio directly 
without using pseudofermions. To begin with, 
we note that the determinant can be written as 
gTr log M g^^^ ^]^g trace can be estimated by the 
Z2 noise which is known to yield the minimum 
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variance Q. Thus the partition in eq. (|^) takes 
the form 

Z = J[dU]e-^''+I '^"^ ^2(iogM)z. 

where Z2 is the Z2 noise for estimating the trace 
of the matrix log M. 

2. Estimating Determinants with Pade-Z2 
Method 

It has been shown that there is an efficient algo- 
rithm to estimate the determinant by expanding 
the logM using the Pade approximation. The 
trace is then estimated by the Z2 noise. This 
is reported by Chris Thron Q. The expectation 
value of the log detM can be written as 

E[logdetM] = E[TrlogM] ^ 1 

where bi and Ci are constants. It is worthwhile 
noting that since we will use this estimate for the 
determinant ratios in the Monte Carlo step for 
the two matrices Mi and M2 whose eigenvalues 
are close to each other, we can do the Pade ex- 
pansion aroud 1 with a small set of bi and As 
long as there is enough memory, the matrix in- 
version with the set of bi and Ci can be done with 
the minimum residual with about 8% overhead in 
CPU time as compared with inversion of the one 
with the smallest diagonal term (the equivalent of 
the quark mass) (^-|7j. Furthermore, Ci are real 
and positive so that the matrix M + Cil is better 
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conditioned then M itself. Finally, we should cau- 
tion that the PZ method does not work if there 
are real negative eigenvalues, because the log has 
a cut on the negative real axis. 

3. Hybrid Monte Carlo with PZ Determi- 
nant 

The algorithm we propose involves two steps: 
1. Molecular Dynamics with Determinant 
Force: 

The molecular dynamics step is retained as in 
HMC, except the pseudofermion force is replaced 
with the determinant force in the evolution equa- 
tions. 



-(y) = -(0)-T[^(0)C/(0)]^ 

U{5t) ^ ea;p[7r(y)(5r]t/(0) 

dlogM St 
+ Z2T[ {St)U{6t)]z2y 



(5) 



Several comments are ready: 

The leapfrog scheme is adopted for this step as 
in HMC to make it reversible and area preserv- 
ing. There is a discretization error iVMi5 0('5T^) = 
0{St^) as in HMC. Like the pseudofermion vari- 
able in HMC, Z2 is kept fixed during the MD tra- 
jectory. 

2. Monte Carlo: 

In lieu of the Metropolis accept/reject criteria, 
we shall us the Kennedy-Kuti linear test |^ 

where U and V are two configurations and R ~ 
g-SH ^ H is molecular dynamics hamiltonian. We 
note that this satisfies detailed balance provided 
R is an unbiased estimate |8| . Since the estimate 
of 5H is unbiased with the Z2 noise, we can use 
the Bhanot-Kennedy algorithm ||] to turn R into 
an unbiased estimate: 



Suppose X is an unbiased estimate. 



(6) 



4. Comparison with HMC with Pseud- 
ofermions 

We shall point out several apparent differ- 
ences between the pseudofermion algorithm and 
the current determinant algorithm and will make 
some preliminary comparison of the two ap- 
proaches. 

I. MD trajectories: 

The pseudofermion approach requires two ma- 
trix inversions — > M'^Mx — x '^'^ each step. 
Whereas, the determinant approach reqires only 
one such inversion in — > Mx = Z2. Accuracy 
of the residual of matrix inversion is needed for 
reversibility for both approaches, i.e. 
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where the average is over all the trajectories. This 
allows one to test the thermal equilibrium and 
make a comparison of the two approaches. 

2. Overhead of Determinant Algorithm in 
Accept/Reject Step: 

The Bhanot-Kennedy unbiased estimate in eq. 
(^), needs on the average e terms in the matrix 
inversion to have independent estimate of SH. It 
may be desirable to have more than one Z2 noise 
in the estimate of the determinant ratio. 

In the Kennedy-Kuti hnear accept /reject test, 
the ratio R should lie between and 2 to be a 
probability. This, however, can be violated in 
practice Q . To make sure that the systematic er- 
rors due to the violation of < i? < 2 is small, say 
less than 0.1%, one needs to have small enough 
step size in the MD trajectories. 

3. Intrinsic Noisiness: 

The ultimate test between different algorithms 
will be conducted by comparing how fast ob- 
servables decorrelate. This will be done next. 
In the meantime, we shall concentrate on the 
issue of intrinsic noisiness. While the pseud- 
ofermion method gives the joint probability of 
gauge fields and pseudofermion fields, the deter- 
minant method has built in it the fluctuation of 
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the Z2 noise. We shall explore the intrinsic noisi- 
ness by computing detM^jdetMi for a fixed pair 
of Ml and M2 separated by heatbath update of 
20 links on a 16^ X 24 lattice with k = 0.154 with 
the two different methods. 

In the pseudofermion method, the determinant 
ratio is estimated by 



detMi/detM2 
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where Xi = -^i^i- In the PZ method , the deter- 
minant ratio is estimated by 



log{ 



detMi , 
detM2 ' 
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We plot in Fig. 1 and Fig. 2 the distribu- 
tion of \og{detMi/detM2) - {\og{detMi/ detKh)) 
for 50 gaussian and Z2 noises respectively. We 
see that the distribution of the pseudofermion 
method is wider than that of the PZ method. 
We also plot in Fig. 3 the Jackknife errors of 
\og{detMi/detM2) due to the two method as a 
function of the number of noise. Again we see 
that the PZ method yields smaller error than that 
of the pseudofermion method. 

Extensive tests are needed to check the effi- 
ciency of this algorithm before we can tell if this 
algorithm will turn out to a useful one. 
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Figure 1. distributions of log{detMi/detM2) — 
{\og{detMi/detM2)) dor gaussian and Z2 noises 
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